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ABSTRACT 


A generalized vaterjet propulsion system in a sub- 
cavitatina hydrofoil craft is considered for design and 
performance estimation. Independent system parameters of 
ee locity ratio, inlet velocity ratio, and nacelle inlet 
diameter to maximum diameter ratio are varied to search for 
the optimum system. The optimum system is defined as the 
minimum total propulsion system weight. The optimization 
scheme utilizes a directed search without the calculation 
of derivatives, and was chosen for its simplicity, versa- 
tility, and rapidity. The ducting system is divided into 
components for head loss calculations. For a given design, 
total head losses in the duct are computed by means of 
experimental data and empirical equations, enabling pum, 
reduction gear, and prime mover desian to be completed. 
Results for a sample craft are included which indicate that 
hydrofoils should be designed about the gas turbine, due to 
the discrete power levels that are available.  Suboptimizing 
the nacelle design may lead to an overall less optimum 
system design. A general lack of information was noted on 
nacelle, cascade corner, and three dimensional diffuser 
design and performance prediction. A FORTRAN computer 
listing and flow charts are included. 
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CHAPTER l - INTRODUCTION 
The use of waterjets for the propulsion of marine 

vehicles at conventional ship speeds has not been consid- 
ered comparable to that of propellers. Although the con- 
cept is quite old, the disadvantages of waterjets have not 
been overcome and the advantages not realized. At the 
speeds considered, the low system efficiency, large system 
weight and relatively sophisticated technology required 
were sufficient to virtually exclude waterjets from com- 
petition with propellers. With the advent of high speed 
craft however, this means of propulsion is receiving a 
renewed interest as a practicable system. Among the pos- 
sible advantages that may be realized are (refs. 1,2): 

a. Reduced noise and vibration 

b. Fewer system components 

c. Elimination of underwater appendages 

d. Elimination of complex transmission machinery 

e. Limited draft applications 

f. Simplification of steering devices 

For hydrofoil craft, the waterjet looks particularly 

attractive. It neatly avoids the problem of conventional 
propellers in power transmission (refs. 3,4). For this 
reason, waterjet propulsion development to-date has been 
primarily concerned with hydrofoil applications. In 


addition, new concepts in ships, such as surface effect 
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ships and captured air bubble craft, are now including con- 
sideration of waterjets as a possible propulsion system. 
With this surge of interest, the need to go beyond the 
practical "Will it work?" engineering design method has 
become paramount. The question has become, "What is the 
m waterjet propulsion system for the craft under consid- 
eration?". Answering this question requires a strictly 
formulated methodology previously not required in make-it- 
work technological levels. No longer is it satisfactory 
Mind some point, any point, above the minimum acceptable 
level. The scope has been enlarged to pinpointing the 
optimum within the entire domain of feasibility. The 
immensity of this task coupled with its relative lack of 
need has hampered the implementation of a methodology. 
Hydrofoils have demonstrated the need, while computers 
provide the invaluable aid for the required computations; 


therefore, an optimization methodology has been developed. 


-11- 





CHAPTER 2 - DESIGN METHODOLOGY 

Родина те thewdesiaDNNWethod, a clear definition 
of the desired objective must be given. In the present 
instance, the definition is the answer to the question, 
"What is meant by an optimal waterjet propulsion system?". 
ШИЕ enoe of an optimum must adequately describe the 
system without bias. It is readily apparent that whatever 
optimum is chosen, it will have a strong effect on the 
methodology employed. The current most commonly used 
descriptor is the propulsive efficiency which is defined 
as the ratio of the power required to the power output of 
the plant (ref. 5). Efficiencies suffer from differences 
m definition and thus vary in their relative importance to 
system optimization. Minimum power required to propel the 
craft at the desired speeds is sometimes utilized as an 
optimization objective with the advantage of being readily 
calculable and understood. A disadvantage is that a gas 
turbine, which is a commonly used prime mover on account of 
its low weight to power ratio, operates most efficiently at 
its maximum power level. Hence, there exists a natural 
trade-off between the most efficient propulsion and the 
margin required for growth in the ensuing years of the 
craft's life. Minimum power calculations are certainly 
necessary but provide a somewhat questionable optimum. In 


an attempt to avoid vagaries of definition and usage and 
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to remove doubts of sufficiency, the goal of optimization 
or the system "cost" may be defined differently. The 
"cost" of any system should reflect the disadvantages of 
that system vithout weighting one disadvantage over 
another. Two disadvantages of waterjet propulsion systems 
mentioned previously are a large system weight, in part 
due to the entrained water, and a low overall system 
efficiency.  Inefficiencies are, however, necessarily 
manifested in increased weight, be it added required 
fuel, another pump stage, or larger ducting. Recall, too, 
that hydrofoils are weight limited vehicles. A unique 
number that reflects these realizations is a suitable 
definition of the weights attributable to the propulsion 
system (ref. 5). The propulsion system weight shall be 
defined as Ele sum Oof all weights of the structural com- 
ponents, prime movers, fuel, pumps, gearbox, ducting, and 
entrained water above the mean waterline. If, for a given 
craft, the displacement is considered constant, then minimi- 
zation of the propulsion system weight is equivalent to 
maximizing the craft payload. Thus, a propulsion system 
shall be defined as optimal if, in addition to fulfilling the 
propulsion requirements, it adds the minimum weight of 
any propulsion system to the vehicle. 

The generalized waterjet propulsion system for 
hydrofoils utilized in this study is depicted in Figure 1. 


It should be noted that this configuration permits 
-13- 
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considerable variation by the changing of a few lengths 
and angles of the ducting. The water enters the nacelle 
under ram pressure and is ducted up the strut and into the 
pump where the energy input is converted to kinetic 
energy in the jet. Thus the thrust reaction is generated, 
providing the craft's propulsion. In equation forn, 

T= 00 (У. сов (В) - Vg) (1) 

In any optimization procedure, the formulation of the 
method requires the knowledge of the parameters to be 
varied. The number and selection of these parameters 
contributes significantly to the applicability of the 
method. The parameters must be necessary and sufficient 
system descriptors, and it is desirable to utilize the 
minimum number required. If M operations are performed in 
the evaluation of one design point and N independent para- 
meters describe the system, then the total number of 
operations required increases as MN. The sufficiency of the 
description of the system is clearly a choice of adequacy. 
The variables that most strongly affect the resulting system 
should be chosen as system parameters while other variables 
that produce only minor variations over a wide range of 
values on the essential design should be defined, or even 
somewhat arbitrarily en, to suit the needs of the design. 

Examination of equation (1) reveals that for a given 
craft and associated mission, only the flow rate and jet 


velocity are unknown. The one-dimensional continuity 
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equation for an incompressible fluid 
AV, = 229; = 0 (2) 

indicates that choice of jet velocity will determine the 
jet area. The first non-dimensionalized system 22 
then, shall be V4/ Vo; defined at the cruise point. Knowing 
the flow rate from the jet, the same reasoning may be applied 
to the inlet which results in the second system parameter of 
Vi/Vg, where Vi is defined in the cruise condition at the 
point of minimum area in the nacelle. The third and 
last system parameter, for reasons that will be explained 
in Chapter 3, shall be D;/Dy, where Dj is the diameter of 
the inlet at the rd о = 

Finally, the minimum amount of information needed to 
specify the craft and its mission must be determined. .A 
typical drag to craft speed curve is shown in Figure 2. 
The "hump" at the lower speeds represents the take-off 
point, i.e., when the hull lifts free of the water and the 
craft is foilborne, and generally occurs at speeds that are 
approximately 50$ of the cruise speed. The sizing of the 
propulsion system then is dependent on both the cruise and 
take-off conditions. The minimum information required is 
the following: 

a. Displacement 

b. Range 

C. Prime mover 


d. Craft Speed, at take-off and cruise 
-16- 
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e. Craft Drag, at take-off and cruise 

f. Depth of submergence of foil 

g. Distance of strut from transom 

h. Distance of pump exit from transom 

i. Height of pump centerline above keel 

Although the prime mover is assumed to be a gas tur- 
bine, knowledge of the particular type of gas turbine is 
important in consideration of both fuel and power re- 
quirements. It is readily seen that the last four input 
requirements define the approximate configuration while 
the first five items describe the craft size and mission. 

Knowing the required ingredients and the desired re- 
sult leaves only the recipe to be set forth. For the 
traverse of an ideal fluid between two points in a closed 
conduit, Bernoulli's equation may be applied such that 

py + loVf + pgz} = p, + lov? + 292. (3) 
where 2 indicates a station downstream of station 1. 
For a real fluid vith losses this becomes 

1272 Е loy? 

Pj * 20V + 092) = р, + 2072 + 092, - Ару, (4) 
where Ар,» = losses between stations 1 and 2. 
Although this equation can be applied only along a stream- 
line and hence cannot be used across the pump, it is 
applicable to the ducting on either side of the pump. The 
continuity equation provides the information that the 
ducting velocities are determined by geometric configura- 


tion. Starting from a known point, e.g. the water surface, 
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and following the water through the duct leaves only the 
losses as unknown. Estimation of the losses may be accom- 
plished based on duct size and shape (ref. 8). This yields 
the needed information for the pump design. 

Following the flow diagram of Figure 3, and recalling 
the previous discussion, the choice of the system parameters 
by the technique described in Appendix A determines the in- 
let and outlet geometry. Knowledge of the jet area determines 
the flow rates at both cruise and take-off. Water is in- 
gested into the nacelle at velocity Vi. The parameter 0, /Ю а 
sets the amount of diffusion possible prior to entrv into the 
strut elbow. Determination of the nacelle lenath provides 
both a drag estimate and the head loss in the nacelle duct 
as described in Appendix B. The strut elbow will exhibit 
large losses due to the low aspect ratio (width/depth) of 
the duct as it conforms to the shape required for the strut 
and due to the relatively large velocities. The use of 
quide vanes will minimize the losses by increasing the 
aspect ratio. Assuming the inlet and outlet areas are the 
same, the losses in the elbow may be calculated by the method 
described in Appendix C. The strut diffuser area ratio is 
set to avoid cavitation in the downstream ducting and to 
reduce the inflow velocity to the pump. Sizing the internal 
flow area permits an estimate of the exterior dimensions 


which are sized up to cavitation limits. For more information, 
-19- 
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see Appendix D. The drags attributable to the strut and spray 
are estimated. The newly calculated nacelle, strut, and 
spray drags are compared to the original estimates. If the 
deviation is excessive, the drag estimate is revised and 

the problem restarted. Otherwise, the head losses in the 
strut due to friction and diffusion are calculated.. Entry 
into the hull from the strut is made via a 90° elbow for 
which the losses are calculated as before. It should be 
noted that this elbow will have a more favorable aspect 

ratio as a result of the strut shape and will require fewer 
guide vanes. Coupled with the decreased velocity, the 

losses will show a decrease from these due to the strut 
elbow. Determination of the ducting required, as described 
in Appendix E, permits estimation of the ducting losses 
therein. Summation of all the losses prior to the pump and 
knowledge of the inflow velocity provide the cavitation 
information in the form of net positive suction head (NPSH) 
to the pump. Adding the nozzle losses permits calculation of 
the pump head required. The method for determining the 
nozzle losses may be found in Appendix G. This is sufficient 
to determine the size and type of pump needed, the details 

of which may be found in Appendix F. The specifications of 
the prime mover are then utilized to specify the gearbox 
parameters of gearbox ratio and power which in turn indi- 


cates the SFC required in the cruise condition. Since the 
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amount of fuel required for the endurance specified is 
calculable, the entire propulsion system is fully defined 
in dimensions, configuration, weights, and particular com- 
ponents.  Summing all the component weights yields the 
desired system weight. This concludes the evaluation of 

a propulsion system weight for the chosen system parameters. 
The optimization routine utilizes the weight information, 
selects a new set of system parameter values and the 

cycle реаіпѕ anew. Ultimately, the minimum weight system 


is located and the process completer. 





CHAPTER 3 - DISCUSSION OF RESULTS 

A sample output is depicted in Figure 6. The first 
page is predominantly a repetition of the input data and 
is printed out upon call to the program. The second and 
third pages are the resulting optimized system data that is 
printed for each configuration of pumps and gas turbines 
considered. Hence the designer can see at a glance the 
pertinent information required for a preliminary design as 
well as compare the alternative configurations that might 
be feasible and competitive designs. 

There are a number of features within the results that 
are not immediately obvious. First and foremost is the 
ھت‎ of results exceeding a system limitation. For 
example, having required cruise shaft horsepower exceeding 
the maximum normal horsepower of the engine. Within the 
program during the optimization, such designs are, of 
course, rejected. For the print out, it was considered that 
knowledge of the resulting "optimum" was desirable, how- 
ever erroneous it may be. One of the forcing influences on 
this decision was the variability of the total drag estima- 
tion and thus the flow rate. If the optimum system re- 
sulted in a value that was at a boundary, a slight change 
in, say, flow rate might cause the value to exceed the 
boondary and hence result in a rejected design. In other 
words the optimum system is not identically reproducible. 


The significance of this fact will be amplified in the 
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conclusions section. This lack of reproducibility results 
primarily from the iterations for drag estimation. Each 
new estimation is compared to the previous estimation and 
accepted if within 5$ of the former value. Otherwise the 
flow rate is adjusted for the new drag and the nacelle and 
strut are redesigned. In the case of the 400 ton craft 
used, 5$ of the total drag is on the order of 4000 pounds 
force in the cruise condition. Thus, starting from the 
same system parameter values but with different initial 
drag estimates may result in a total drag differing by 5$ 
in the output with corresponding differences in the system. 

Secondly, the resulting design is obviously constrained 
by the ducting configuration used. Hence direct comparison 
with existing craft is not valid until the configuration 
is modified to conform to that of the craft. 

At the inception of this project, it was believed 
that the duct area ratio of the strut was a system parameter 
to be varied independently. Clearly a lower bound on this 
parameter is that amount of diffusion required to avoid 
cavitation in the hull ducting. It was felt, however, 
that the pump, gear and fuel weight calculations might 
reflect a lower overall weight if the pump inflow rate was 
reduced. The results did not justify this assumption. The 
system unerringly chose the strut with the minimum required 
diffusion, which meant minimum.weight. The principal 


reason for this result is that the pump design is in terms 
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of total head. Hence the pump effectively never "sees" the 
amount of static head presented to it. The pump is affected 
only by the differences in ducting losses which are re- 
flected in both the NPSH and pump head. Thus the net effect 
is of secondary importance and is not significant enough 

to change the area ratio of the strut. For this reason the 
area ratio of the strut ducting was rejected as a system 
parameter and is presently determined as required by 
Cavitation limitations or pumo inlet area. 

The initial nacelle design method chose the largest 
inlet to maximum diameter ratio permissible by exterior 
Cavitation. This ratio was tempered by the amount of 
diffusion required to avoid cavitation in the strut elbow. 
Assuming the same inlet velocity ratio, would more diffusion 
than the minimum required permit an-overall lesser weight 
system? In other words, does the optimum system require 
the optimum nacelle or might an off design nacelle be 
used resulting in a lesser weight system? The use of the 
nacelle diameter ratio, DIDM, as a system parameter pro- 
vides the opportunity to investigate this question. The 
current results are somewhat inconclusive due to the influence 
of other factors.  Tentatively, it is suggested that for some 
designs, the use of DIDM as a system parameter is justi- 
fied and results in a better system design. It is 
recommended that investigation of this question be pursued 


for determination of what constitutes the basis of nacelle 
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design in regards to system performance. 

Perhaps the most significant conclusion is that the 
choice of engine is extremely important in the resulting 
optimum system. The first reason for this is that the 
gas turbine is presently available only at discrete power 
levels. Since fuel weight is large relative to the other 
system veights, the operation of a gas turbine at a point 
significantly off its design point will result in a large 
increase in fuel weight and thus system weight. The net 
effect is that, given the opportunity to avoid the penalty 
in fuel weight, the system tends towards the maximum SHP 
permitted by the engine, and will drive the jet velocity 
ratio as hiah as possible to achieve this SHP. This also 
results in a low flow rate, low drag and low propulsive 
coefficient. If the engine is slightly large for the craft, 
then the results vill show a high jet velocity ratio, SHP 
per engine at the maximum permissible level, which may be 
limited by take-off or cruise SHP, and a low propulsive 
coefficient. If the engine is so large that there is no 
power limitation, the jet velocity ratio will be very 
high and the SHP per engine relatively low resulting in a 
very low propulsive coefficient due to large total SHP 
and a high system weight. The effect of an engine that is 
barely adequate is a lowered jet velocity ratio in an effort 
to reduce the SHP requirements. The trade-off here is the 


large flow rate required to compensate for the low jet 
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velocity ratio and the large drag that results from 
attempting to accomodate the flow rate in the nacelle and 
strut. Thus SHP is being forced up and the system may never 
achieve an acceptable horsepower level. 

This power level sensitivity is so strong that it is 
abundantly clear that hydrofoils should be designed 
around a particular gas turbine. Choice of a desired 
displacement and the use of an adequate prime mover may 
result in a poor design and inefficient use of the gas 
turbine. Choice of the gas turbine and permitting the 
displacement to be chosen at the point of best engine 
operation vill result in a more efficient and more effective 
design. Thus it is recommended that another version of this 
program be created that chooses the best displacement 
given the engine. It is also suggested that perhaps the 
concept of hydrofoil craft design be revised to accomodate 
the fact that the engine choice determines the best design. 
Figure 7 shows the minimum weight ratio for one LM2500 
and one pump with a lift to drag ratio of 11.3. The 
results show the optimum weight ratios with and without 
fuel do nee coincide. The propulsive coefficient is drawn 
to show the trade-off between maximizing SHP and maximizing 
the propulsive coefficient. 

The choice of an engine also affects the reproducibility 
of the optimum system design. If the engine is smaller 


than desirable, VJVO is driven as high as the SHP permits. 
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When the optimum system parameters have been found and the 
optimum system is to be regenerated, the design process is 
reinstituted. During this final design, iterations in 
Шагы агас are again required to determine the final flow 
rate and size of the system. Thus the net result may be 

a system slightly different from the one originally found 
with least weight. And in the case of a system already on 
a boundary, the resultinq system may exceed the limits 
placed on the optimization process. It is to be expected 
that the deviation will be small but nonetheless present. 
The results in such a case should be interpreted as optimal 
at the boundary and not an unsatisfactory design. If the 
deviation is significant and the other system values 
reflect it, the interpretation should be that the system 
was unable to be optimized within the limits. The most 
common example of this occurs when the cruise or take-off 
SHP of the engine is insufficient to meet the craft re- 
quirements. 

The effect of range on the system design is shown in 
Figure 4. The LM 2500 is a power limited engine for the 
four hundred ton craft at forty-two knots. The FT4C-2, 
however, is so large that maximum power is never attained in 
this craft. The effect of lower SFC and operation at 
maximum power favors the LM 2500 at the longer ranges. As 
it should be, the weight ratio without fuel is essentially 


constant for both engines, although the weight ratio for the 
Dip 





LM 2500 is higher, presumably on account of the power limita- 
pron. Similarly, the jet velocity ratio remains constant 

as does the flow rate, but the inlet velocity ratio tends 

to increase with increasing range. The nacelle diameter 
ratio, on the other hand, shows no consistent pattern. 

Te effect of variation of displacement by 10% is shown in 
Figure 5 for the FT4C-2 for one and two pumps. 

The results show a surprising sensitivity to the 
starting values of the system parameters. When arbitrarily 
deviated by twenty per cent from the present values, the 
results were unpredictable. The current values, however, 
have shown good stability with test craft ranging from 
fifty to seven hundred fifty tons and from forty to fifty 
knots. The reason for this apparent instability is un- 
known. The step sizes appear to be acceptable at any size 
above a threshold value. Some conditions vill not permit 
location of the optimum if the step size is too small. The 
current step sizes, although generally too large, have per- 
formed well for the test craft above. The probable rationale 
is that the large step size permits a search over a wider, 
mere extreme range of system values, but in the event that 
the step is too large, the step size will rapidly be 
reduced to an acceptable level. | 

The lack of experimental data and analytical techniques 
has been noted in several areas. The nacelle design is a 


case in point. Although nacelles have been studied exten- 
— 





sively by the aeronautics industry, the application of 
nacelles to a water environment has been neglected and 
consequentlv, design considerations have not been defined. 
The design of cascaded corners has somewhat more experimental 
data for loss estimation, but the scatter of the data and 

the lack of definitive design information render elbow 

design a vague approximation at best. Three dimensional 
diffusers are likewise an unexplored domain. The effect of 
changes of ductina shape on losses is not treated in the 
literature for straight pipe. Loss estimation for changing 
shape in elbows is treated, but only for rectangular forms. 
It is apparent that the prodram results contain a great deal 
of information, but it is also evident that the interpretation 
of.the results can be improved by a more thorough study of 
runs for different craft, configurations, and missions. А 
study of the effects and interactions of such changes will 


provide a finer understanding of the import of the results. 
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CHAPTER 4 - CONCLUSIONS AND RECOMMENDATIONS 

A method of optimizing a waterjet propulsion system 
for hydrofoils has been developed and tested vith a four 
hundred ton hydrofoil at a cruise speed of forty two knots. 
Other craft have verified the importance of the selection 
of the gas turbine for a particular craft. It is 
recommended that a version of this program be adapted to 
determine the optimum displacement for a given engine. 

Although originally considered important, the amount 
of diffusion in the strut should be the minimum required to 
avoid cavitation in the ducting or the minimum required by 
the pump, whichever is greater. The desian of the nacelle 
need not be at the limits of cavitation, but should be 
integrated into the total system design to provide a lesser 
weight propulsion system. It is recommended that a study of 
nacelle design be initiated to determine the interaction 
of the nacelle and the propulsion system. 

The initial values of the system parameters may affect 
the operation of the optimization procedure. A determination 
of the cause of this lack of stability should be made. The 
magnitude of the step sizes of the parameters has little 
d rect on the program provided they are sufficiently large. 

A general lack of data and design information is 
prevalent for fluid flow in closed conduits. This is 
especially noted for ram inlet desiaqn, cascaded corners, 


three dimensional diffusers, and changes of pipe shape. 
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The program provides a unique facility for the 
preliminary designer to rapidly ascertain the optimum 
waterjet propulsion svstem design and performance. Many 
results should be studied to deduce the complete real 


world interpretation of a computer generated design. 
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APPENDIX A - OPTIMISATION 

The selection of an optimization procedure represents 
a key consideration in the overall design method. For one 
thing it may require computations to be made solely for its 
usage, as in gradient methods. Then, too, it may require 
a cory program to tell it what to do, as in random 
search techniques, or an adaptive direction capability may 
be self-contained (ref. 18). It is not the intent of this 
investigation to find or invent the "optimum" optimization 
procedure, but the significance of its usage should not be 
minimized. 

The general consideration of how to optimize has a 
wide variety of tools to aid in the achievement of the 
answer. One of the more common methods utilized is some 
variant of the gradient search (ref. 19). The rationale 
is an intuitive one, based on the notion that if the goal 
is the attainment of a maximum, for example, the quickest 
way to go is in the direction of the steepest slope. This 
is not unlike a blind man's method of climbing a mountain. 
Like the blind man, the gradient search can encounter 
difficulties by following this philosophy. The assumption 
of a finite slope in the area of interest not only requires 
the absence of steps on the hyversurface, but imposes the 
more stringent limitation that the slope be calculable and 
containable within the computer. This somewhat limits the 


types of systems that can be optimized by this technique and 
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and exacts a stiff penalty for some numerical inaccuracies. 
Butter, assuming that the system to be optimized permits 
slope calculations, it is not known a priori that such 
calculations provide sufficient information to compensate 
for the time spent in making them. In the present system 
Шетел, discrete steps are virtually guaranteed due to the 
different pumps and gearboxes used. Since computation of 
the propulsion system weight is not done by means of an 
analytic function, slope calculations become as lengthy as 
the weight calculations themselves. For these reasons, 
gradient search techniques are rejected as a possible 
optimization tool for this design problem (ref. 20). 

| On the opposite end of the optimization spectrum lies 
the random search. The principal advantages of a random 
search are the ease with which it may be understood, hence 
used and programmed, and the simplicity of its computation. 
No slopes or step sizes need to be calculated or estimated, 
merely the values of the function at the points of interest. 
This is also its major disadvantage. 

The random search also has no greater likelihood of 
finding a local optimum than it does of finding the global, 
and, hence desired optimum. Once the capability of 
"directionability is gained, local optima are just as likely 
to be found as the global optimum as is the actual optimum. 
Ihis is true of all directed searches. The only way to 


avoid this unpleasant result is to use random starting points 
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or to use a systematic search over the entire hypersurface. 
A unimodal surface, of course, does not present this 
problem. 

Since the random search utilizes no other additional 
information than the function value, it is highly in- 
efficient with the result of greater time spent in the 
search. All optimization is done in comparison to the 
previous best. Thus the simplicity of its procedure is 
gained only at the sacrifice of the large number of evaluations 
that must be made to assure a reasonablv high certitude of 
attaining the optimum. This becomes especially significant 
when the time of calculation of a single function value is 
lengthy. In the design and performance estimation of a 
propulsion system, a lengthy calculation is likely to be the 
case, in part due to the iterative loops required (refs. 21, 
22), 

The solution is clearly somewhere in the middle of 
these two approaches, utlizinc the best features of each and 
minimizing their disadvantages. It should have the ability 
to discern direction without the requirement of additional 
calculations. One such method is the pattern search utilized 
in this design (ref. 23). The basis of this method is akin 
to that of the gradient search but is modified in its ap” 
proach. Again based on intuition, the pattern search depends 
on the supposition that the search effort should be expended 


in the direction of improvement. Note that this is a 
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qualitative judament rather than a quantitative one, as with 
gradient methods. The criterion for a sense of direction 
lies in the success of an evaluation, not in its measure 

of success, which in itself would be somewhat arbitrary. 
Thus, while gaining adaptive information about the location 
of an optimum, the disadvantages of the gradient search are 
neatly bypassed and the random search is effectively uti- 
lized in a more efficient manner. 

Following the flow diagram of Figure 8, the pattern 
search uses the starting point as the initial basepoint. 
From the basepoint exploratory moves are made. That is, 
each coordinate is varied from the basepoint position the 
amount of the step size and the function evaluated at that 
point. Each exploratory move is examined for possible 
improvement in the function value. If improvement cannot 
be made in any direction, the step size is decremented and 
the exploratory moves resumed from the basepoint. If a 
position is located with a better function value than the 
value at the basepoint, the second type of move, a pattern 
move, is made. A pattern move is a simultaneous change in 
all dimensions of length equal to the distance from the 
current basepoint to the present position. If improvement 
is made: in either the pattern move or the sequential explor- 
atory moves, another pattern move is made. Lack of improve- 
ment resets the basepoint at the position of the best 


function value and enters the cycle again with exploratory 
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moves from the current basepoint. 

This ability to seek, find, and establish a pattern 
gives this search an ability to move rapidly in the best 
direction. It also has the capability of shifting direc- 
tion or modifying the pattern when required. Step changes 
in the hypersurface present little problem since slopes are 
not utilized. With these features the pattern search can 
move quite readily along steep ridges or narrow valleys 
which are common "Waterloos" of other search methods, 
especially aradient methods. 

Perhaps it should be mentioned here that the intent 
is not to totally degrade the gradient search which has its 
own attributes and an important usefulness in its own right. 
It is not without shortcomings, however, that might present 
a serious detriment to the user who is looking for a 
generalized search method. 

The original form of the pattern search as proposed in 
reference 23 has been modified to conform to the needs of 
this design problem. An additional modification has been 
made that is considered an improvement on the initial 
scheme. During the exploratory move phase of the search, 
each dimension was initially incremented by a step size and 
the function evaluated. If the move was a success, the point 
and its associated value were saved and the next dimension 
explored. Lack of success meant decrementing the position 


by a step size and evaluating the function. Again a success- 
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ful move was followed by retaining the position and value and 
exploring the next dimension. An unsuccessful move was 
ignored and the next dimension NEYE na HEE. The modification 
implemented is to vary the dimension first in the direction 
of the last improvement rather than simply incrementing it. 
Thus anticipating success, time would be saved by not 
checking a poor direction. The resulting advantage of this 
modification is, of course, very dependent on the function 
being evaluated, but is in keeping with the philosophy of 


the search and is accomplished at a very small cost of time. 
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APPENDIX B - NACELLE 

The details of this section may be found in reference 
9, from which the following Comments are drawn. 

The selection of the system parameter values sets the 
ratio of the inlet diameter to the maximum external diameter, 
DIDM. This effectively limits that maximum amount of diffu- 
sion that may be carried out in the nacelle. The present 
design assumes that the maximum amount of diffusion allowed 
in the nacelle duct is utilized. The nacelle is first 
checked to ensure that it does not cavitate externally. 
Prediction of external cavitation rejects that selection of 
system parameter values. If cavitation is not predicted, 
Bue inlet is checked for cavitation at the point of maximum 
flow, which generally occurs in the take-off condition. If 
cavitation was predicted internally, auxiliary inlets are 
sized to reduce the inflow rate to just below the point of 
cavitation. The length of the nacelle is then suboptimized 
to be the length at which the total power loss due to drag 
and diffusion in the duct is minimized. The total head 
loss, made up of the lip loss and the diffuser losses may then 
be calculated. Likewise, the drag coefficient may be 
calculated from Hoerner's approximation for stramlined 


bodies (ref. 10), based on wetted surface, 


3 
CanCe, (141.5 (2) 747 (Br) *) 
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APPENDIX С = ELBOWS 

All elbows in the ducting are treated in identical 
fashion, although individually each elbow may be very 
different. The inlet and outlet areas for any bend are 
assumed to be the same. An elbow is considered fully defined 
by the knowledge of four quantities: width, depth, cross- 
sectional area, and radius ratio. The radius ratio, Е 
is the ratio of the radius of the centerline of the bend to 
the duct radius, as depicted in Figure 9. Since the flow 
rate, Q, has been previously determined, the first three 
items are sufficient to calculate the aspect ratio, shape, 
and the mean velocity of the fluid in the elbow. 


For Re (£9)? > 91, Ito's expression for the head loss 


coefficient in a circular bend is (ref. 11): 


K,2.00241a0Re7?* 1 * (E. ) 9. 9* (C. 1) 


О 


К 
a is a numerical factor dependent on 9 and Ens" For the 


R 
раъд о = 
сасе О Ө = 90? and Бос 19.7 
гс 
a=0.95+17. 2 (Zo) . (С.2) 
The loss coefficient then is 
В: ı-1,96 -0.17,R 10,84 
Ke 0.00241010.95+17.2 (50) ] Re (во) 


For а given duct and considering 0, V, and R constant, 


partial differentiation of Kr with respect to RO leads to 
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= = 0 = (Fo) opt но 
Since к. E = + CRO for these conditions, this result 
verifies the intuitive feeling that a minimum exists due to 
a trade-off between friction and bend losses as the number 
of guide vanes in the elbow varies, as depicted in Figure 10. 

If the given radius ratio is greater than one, the 
losses in the elbow may be reduced by the use of splitters 
which aid in the prevention of separation. In this context, 
splitters are defined as guide vanes of unequal length, 
unequal spacing and possibly unequal curvature (ref. 12). 
Turning vanes, on the other hand, are always equal in length 
and curvature, although the spacing between vanes is not 
necessarily equal. If the splitters are considered in order 
to subdivide the elbow into smaller elbows, the minimum 
loss will occur when the ratio of the inside radius to the 
outside radius of any of the subdivided elbows is constant. 
Referring to Figure 6, this mav be written as 

(сі) = (Бір 
YO A 

where n is the number of subdivided elbows, i.e. the number 
of splitters + 1. Figure ll depicts a duct with one 
splitter. Clearly, 


R = RO 


қ 
И 


у = В+ RO 
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Hence, the desired radius ratio, 4.3, may be substituted for 


(хода) and the resulting equation Solved for n, or, 





ENS “ia 
Ro 
Ln 
R R 
ILE Ser 
n= RO н = 2.12 Ln IRO 1 (СТА) 





37 


Ln Беті 


The number of splitters determines the spacing. The head 
loss coefficient may be calculated for each subdivided 
elbow and the results averaged over the area which will 
yield the average head loss coefficient for the elbow. 


For R < l, equation C.4 has no meaning. For this case 


RO 
the total head loss coefficient was computed as the sum of 


the friction loss coefficient and the loss coefficient from 
Figure 12 which is an empirical graph tor thin, circular arc 
turning vanes from reference 12. 

In either case the results are valid up to a Reynolds 
number below that expected in this design. The head loss 
coefficient includina friction is then corrected for the 


higher Reynolds number by (ref. 13) 


Ea 


where by curve fit, 
x = 1.0057-0.16892Ln(Re x 10 °) +0.0145885[Ln(Re x 10”°)]? 


-0.0016948(Іл(Ре х 1075)13 (6555) 
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APPENDIX D - STRUT 
The head losses in the strut diffuser are caused by 
friction and the area expansion ch thickens the boundary 
laver. The friction coefficient 15 computed in the same 
manner as for straight pipe. The expansion loss coefficient 
is 


ہم 
K, = C(1- AS (D. 1)‏ 


C is dependent on the diffuser equivalent angle, 286, and 

is interpolated from the graph in Figure 13, which is for 

a two-dimensional, straicht-valled diffuser. For the flow 
area, the thickness to chord ratio of the strut is held 
constant while the area is increased. Hence the assumption 
Mewtwo dimensional diffusion is not strictly valid but is 
considered a conservative estimate. 

The strut dissipates energy externally by a drag term 
attributable to the body and a term resulting from the spray. 
The drag coefficient from Hoerner (ref. 10) is based on the 
planform area and is approximated by 

Ca. = 2Cg, [1*2 (t/c) *60 (t/c) *] (D.2) 

The spray drag coefficient, likewise based on planform area 
15 

я 20/5: 06/46) (рез) 

This approximation is valid for Froude number, based on the 


chord length, greater than three (ref. 14). 
sip 





APPENDIX E 
Pump Inlet Piping 
The ducting from the hull elbow to the pump inlet is 
a straight pipe, a 90? elbov and a transition piece. The 
length of the straight pipe depends on the number of prime 
movers. In equation form, 


Length = Beam 
Number of prime movers + 1 


The transition piece is shaped according to the number of 
pumps under consideration. Since there are two struts 
ducting water into the hull, assuming a canard foil 
arrangement, and only one, two, or four pumps are consi- 
dered, the transition piece will be a junction, straight 
pipe, or diverging branch. 

All straight pipe loss calculations are made usina the 
standard Darcy-Weisbach equation 


m= f y? (E. 1) 


The Moody friction factor, f, is computed using the closed 
form approximation of Techo et al (ref. 15), 

f = [.86859 Ln(Re/1.964Ln(Re) - 3.8215))172 (E.2) 
Ro is based on the mean hydraulic diameter which is defined 
as 


ec Mp (E. 3) 


This approximation for the friction factor is accurate to 


within .1$ over the range of Reynolds number from 104 «о 
CE 





Пони 108. 


For both the junction and divergence, the average 
velocity is assumed to be the same before and after the 
transition piece and the piping is symmetrical about the 
single pipe as shown in Figure 14. This requires that the 
smaller piping cross-sectional area is 1/2 that of the 
larger piping. 

The mixing loss coefficient in the junction, without 
friction loss, is computed from 

K 4 = ] + Х ~ 2со5(Е(ф)) (Е.4) 
where f($) з 1.4ф - .005836? 
апа ф 15 the angle the small branch makes with the larger 
duct (ref. 16). In the cases examined for this design, 6$ 
was set to zero resulting in only friction losses in the 
Пе зоп. 

The divergence loss factor A was computed from Figure 
15 which is adapted from reference 16. The angle of the 
divergence, of course, depends on the pump length. It 
should be noted that the divergence configuration requires 
another elbow, equal in angle to that of the divergence, 
to enter the pump. 

In all three transition pieces the size was set as 
required, depending on the pump le". If the pump's 
length requires a zero or negative transition piece 
length, the transition piece is neglected and the losses 


computed only for the athwartships length and the pump elbow 
| =58-— 








turns the ducting aft. 
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APPENDIX F 
Pump 

The details of this appendix may be found in reference 
17, from which the following comments are dravn. 

There are two basic types of pumps Bu s in this 
design: а multistage, double impeller centrifugal design 
and an axial pump with an inducer impeller. The variations 
occur in the number of stages used. The pump design 
utilizes the flow rate and Thoma's criterion is defined as 


о = NPSH (Eel) 
Pump Head 


where the net positive suction head, 


NPSH = V6 + ha - h, -hegjey - Ah, (F.2) 
га 
апа 
Pump Head = у2 - у2 + helev + Ah, око 
2g 


The axial pump's first stage is an inducer impeller 
with up to two additional stages available if the amount 
of head is sufficiently high. Assumption of a specific 
Speed yields the desian pump speed and approximation 
to the off-desian pump speed, which in this case is the 
Sauce condition since the pump is basically designed 
for take-off conditions. The RPM of the prime mover deter- 
mines the gear ratio and calculation of the shaft horse- 
power permits an estimation of the gear weight. Finally, 


a fuel weight may be calculated for the range specified in 
E 0s 





Ehe cruise condition. 

The centrifugal pump is designed in the same way except 
that the number of double suction impellers is increased, 
thus increasing the RPM, until the combined weights of the 
pump, dry and wet, gearbox and fuel are at a minimum. 

The two resulting pump designs are then compared on 
the total pump, fuel, and gearbox weights and the least 


weiaght system 1s selected for the propulsion plant. 
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APPENDIX G 
Nozzle 
For the purpose of calculation the nozzle vas consid- 
ered a straight-walled nozzle. Assuming the losses are of 
the .form 


h = Є г у? (G. 1) 
n = = == 
р 
and noting the friction factor is approximately constant 


and the velocity essentially uniform, this may be written 


Ah = f AL V? 
> р 24 
Or 
T 
hn = f fil dL 
0 2а 
Substituting D - р. + Ltan (y) 
V= 40 
"D? 
results in 
16 2 
na = £ 24/2 [lp] J (6.2) 
De7D4 DT 2а 


which is used for the head loss estimation in the nozzle 
when р, » D. SL, (rel 2). When < 1, the head loss 
is assumed to be 1.5$ of the dynamic head. 

In the event that the distance of the pump exit from 
the transom is greater than 2 times the height of the pump 


centerline above the keel, the nozzle is assumed to exit 


through the underside of the craft. The latter distance 
=62=— 





is that if the nozzle were directed along an axis at 45? 


from the 


vertical. 


In every instance, however, the nozzle is assumed to 


be depressed at the optimum angle to reduce the lift 


required which would allow the use of smaller foils and 


hence reduce drag. As before, 


| = 


Rewritten this is also 


5 = 


where 


T 
M 


Equating 


о0 (У.сов (8) - u (1) 
(= т) е (С. 3) 


= pQV.sin (8) 


equations (1) and (G.3), and solving for Е gives 


А = (pQV, + WD) [cos (8) + Dsin (8)]71 
L 


L 


The minimum energy input to the jet occurs when V; is at a 


minimum, 


oV. 
387 


which is 


= 0 => (pQV, + WD) (рсов(8) - 51п(8)) 
L L = 0 
(cos(8) * Dsin(8))" 
L 
=> р т TAN | D ] 
É (G.4) 
the optimum angle for nozzle depression (refs. 4,6). 
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USING ND DERI vATI VES. 
REF: HOKE Є 
JEEVE Se JACM, APRIL 
1961 +PP 212-229. 
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Label 


PARMS 


DRAG 


FLOW 


LIST OF “SYMBOLS 


USED IN COMMON 


Variable 


Name Subscript 


VJ VO 
VIVO 


DIDM 


TDRAG (I) 1 


STRTD (I) 


POD (I) 


SP RAY (I) 


REST (I) 


VO (I) 


ТВІМ (1) 


О (1) 
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Description/Usage 

System parameters 

Jet velocity ratio 

Inlet velocity ratio 
Inlet diameter to 
maximum nacelle diameter 
ratio 


Drag estimates 


Total craft drag at 
condition I, pound force 


Cruise condi ti on 
Take-off condition 
Evaluation conditions 


Strut drag at condition 
I, pound force 


Nacelle drag at condition 
I, pound force 


Spray drag at condition 
I, pound force 


TDRAG (I) -SPRAY (I) -STRTD(I) 
-POD(I), pound force 


Craft velocity at condition 
I, foot per second 


Craft trim angle at con- 
dition I, degree 


Flow characteristics 
Volume flow rate at con- 


dition I, cubic foot per 
second 





ELBW 


CHARS 


Variable 


Мале Subscript 
AIN 
AJET 
AREA(ICOMP) 
VJ (1) 
VI (I) 
XK ( IELB) IELB 
l 
2 
3 
4 
RO (IELB) 
THATA (IELB) 
wIDTH 
DEPTH 


TYPE (3, TYPE) 


WGTS (LS,LC) 


= Ре 


Deseription/Usage 


Inlet area, square foot 


Jet area, square foot 


Ducting cross-sectional 


area at station ICOMP, 
square foot 


Jet velocity at condition 
I, foot per second 


Inlet velocity at con- 
dition I, foot ver second 


Elbow data 


Radius ratio of elbow 
LELB 


Strut elbow 

Hull elbow 

Pumo elbow 
Divergence elbow 


Duct radius at elbow 
TEEB, foot 


Angle of bend of elbow 
IELB, degree 


Width of duct at elbow 
inlet 


Depth of duct at elbow 
inlet 


Contains name of elbow 
shape ITYPE 


System characteristics 


Weight of component LC, 
pound 





"m 
0 














mu 


Variable 
Name 


CGS (LG,LC) 


Subscript 
LS 


10 
11 
12 
ЕЗ 


14 


15 


LG 


2952 


Зате‏ / تہ ۶1 ۶1ء 5 6ق 


Structure . 


Water (or fuel) 


Nacelle 

Strut elbow 
Strut diffuser 
Hull elbow 
Athwartships length 
Pump elbow 
Transition piece 
Pump 

Nozzle 

Reduction gear 
Fuel 

Prime mover 


Lift from nozzle 
depression 


Total weight 

Spare location 

Centers of gravity of 
component LC, excluding 


fuel, gearbox, or prime 
mover, foot 


Structure vertical center 
gravity (from keel) 





Label 


H20 


Variable 
Name 


DELH (I,ICOMP) 


СС5Х 


CGS2 


TEMP 


РМ 


RHOW 


Subscript 


2 


ICOMP 


105415 


-пе- 


Description/Usage 


Structure longitudinal 
center of gravity 
(from transom) 


Water vertical center of 
gravity, (from keel) 


Water longitudinal center 
of gravity, (from transom) 


Total head loss up to and 
including component ICOMP 
abecondition I, foot 


Same as defined except 
Nozzle head loss only 


Total system head loss 
(including elevation), 
DELH(I,9)=DELH(I,8)+ 
DELH (1,7) 


Spare locations 


Total longitudinal center 
of gravity of system 
excluding gearbox, prime 
mover and fuel, from 
transom, foot 


Total vertical center of 
gravity of system exclud- 
ing gearbox, prime 

mover and fuel, from keel, 
foot 


Sea water (3.5% salinity) 
properties 


Temperature, degree 
Fahrenheit 


Vapor head, foot 
Density, pound force 


second squared per foot 
to the fourth 





Label 


TOLER 


РБОВ 


Variable 
Мапе 


GNU 


HA 


DELTA 


GERAT (NSTG) 


SHP (I,NSTG) 


RPM(I,NSTG) 


PERF (L,IENGN) 


ETAP (I,NSTG) 


subseript 


- E 


pescription/Usage 


Viscosity, foot 
Squared per second 


Atmospheric head, foot 


Check for optimum system 
Pump system 


Gear ratio required for 
pump NSTG 


Shaft horsepower re- 
quired at condition I for 
pump NSTG, horsepower 


Axial speed of pump NSTG 
at condition I, revo- 
lution per minute 


Prime mover IENGN 
Characteristics (L) 


Maximum normal horsepower 
at design speed, horse- 
power 


Maximum intermittent 
horsepower at design 
speed, horsepower 


Specific fuel consump- 
tion (SFC) at desian speed 
and maximum normal horse- 
power, pound fuel per 
horsepower hour 


Design speed, RPM 


Prime mover weight, with- 
out auxiliaries, pound 


Efficiency of pump NSTG 
at condition I 








SHEP 


NACELLE 


Variable 
Name 


DESP 
RANGE 
BEAM 


HS 


HE 


HCL 


XLS 


XLPE 


X LP 


DRAT 


DM 


AI 


AIAUX 


ELEXT 


ELENT 
oe 


Subscript 


Description/Usage 

Same as defined except 
Reduction gear 

Spare locations 

Craft characteristics 
Craft displacement, pound 
Endurance, nautical mile 
Ерос 


Beam, 


Depth of submergence of 
foil, foot 


Height of pump centerline 
above keel, foot 


Height of pump centerline 
above mean water, foot 


Distance of centerline 
Of strut root from tran- 
som, foot 


Distance of pump exit 
from transom, foot 


Length of pump, foot 
Nacelle characteristics 
Diameter ratio, DI/DM 


Maximum external diameter, 
foot 


Inlet area per nacelle, 
square foot 


Auxiliary inlet area per 
nacelle, square foot 
Length of forebody, foot 


Lenath of lip, foot 





Label 


CONST 


ВЕС 


INDEX 


Variable 
Name 


ELAUX 


ELDIF 


ELN 


PI 


RHOD 


те 


PL 
el 


CPM 


IEVAL 


Subscript 


«0 


<0 
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Description/Usage 


Length due to auxiliary 
inlet, foot 


Length of diffuser, foot 
Lenqth of nacelle, foot 
Constants 

3.14159265 

Acceleration of aravity, 
32.174,foot per second 


squared 


Density of steel, 480, 
pound per cubic foot 


, Strut characteristics 


Thickness to chord ratio 
Thickness at root, foot 


Chord at root, foot 


. Thickness at tip, foot 


Chord at tip, foot 


Chord at flying water- 
line, foot 


Indices for program 
control 


Design/evaluation index 


Evaluate at |IEVAL| points, 
no design 


Design at cruise and take- 
off conditions only 


Design at cruise and take- 
off conditions and evaluate 
at |IEVAL| points 





Label 


Variable 
Name 


IEOPT 


ESTERI 


NUAB 


ТЕМОМ 


Subscript 


>0 


2 


2+ | ТЕМАЦ | 


10 
11 


12 
=@0- 


Description/Usaae 


No equipments or config- 
uration entered; total 
system desian/evaluation 


Specific equipments or 
configuration entered 
(not yet implemented into 
program) 

Initial proaram condition 


Design and/or evaluate 


Evaluate only 


. Final program condition 


Design only 
Design and/or evaluate 


Prime mover type 


| TF 35 


TF 40 

Proteus, 1500 rpm 
Proteus, 1000 rom 
Тупе ЈА 

Tyne ІС 

ЕТ12А 

LM 1500 

LM 2500 

FT4A-2C 

FT4A-12 


ЕТ4С-2 





Label 


ITABL 


CDRAG 


WEGT 


Variable 


Name Subserint 

ITYPE 
1 
2 
3 
4 

ICOMP 
l 
2 
3 
4 
5 
6 
- 
8 
9 

NPUMP 

NGT 

CSTRT(I) 

CPOD(I) 

СБРЕУ (I) 
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Description/Usaae 
Elbow shape 

Ellipse 

Circle 

Rectangle 

Square 

Component index 
Nacelle outlet 

Strut elbow 

Strut diffuser exit 
Hull elbow 
Athwartshins length 
Pump elbow 

Transition piece 

Pump inlet 

Nozzle throat 

Number of pumps 

Number of das turbines 
Intervolation parameters 
Computed drags 


Computed strut drag at 
condition I, pound force 


Computed nacelle drag at 
condition I, pound force 


Computed spray drag at 
condition I, pound force 


Pump, gear and fuel weights 


| 








0 
жи 


د 


. нэ 








Label 


PUMM 


Variable 
Name 

XWD (NSTG) 

XWW(NSTG) 


XWG (NSTG) 


XWF (NSTG) 


QQ (I) 


D1S (NSTG) 


D2S (NSTG) 


XNS (NSTG) 


SM(NSTG) 


PLP (NSTG) 


NSTG 


Subscript 


202- 


Description/Usage 


. Pump NSTG dry weight, 


pound 


Pump NSTG wet weight, 
pound 


Gearbox weight associated 
with pump NSTG, pound 


Fuel weight associated 
with pump NSTG, pound 


Pump characteristics 
Flow rate per pump at 
condition I, cubic foot 


per second 


Inlet tip diameter of 
pump NSTG, foot 


Exit tip diameter of pump 
NSTG, foot 


Specific speed of pump 
NSTG, cfs units | 


Suction specific speed 
of pump NSTG, cfs units 


Length of pump NSTG, foot 
Indicator of pump type 


Axial pump with single 
inducer impeller 


Axial pump with inducer 
and one impeller stage 


Axial pump with inducer 
and two impeller stages 


Centrifugal pump with max- 
imum of ten parallel 
double suction impellers 


Spare location 





Variable 


Label Name Subseript 


SHI (NSTG) 


AI 


HEAD 


HPP (I) 


HSV (I) 


THOM (I) 


PHI (NSTG) 


WF 


WG 


RARA 


IFUEL 


e TRUE. 


FALSE. 
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Description/Usaae 


Head coefficient of pump 
NSTG 


Number of parallel double 
suction impellers for 
centrifugal pump 

Number of impellers, not 
including inducer, for 
centri fugal pump 

Pump flow characteristics 


Pump head at condition I, 
foot 


Net positive suction head 
at condition I, foot 


Thoma's cavitation index 
at condition I 


Flow coefficient of pump 
NSTG 


Working variable for fuel 
weight, pound 


Working variable for gear 
weight, pound 


Unlabeled common variables 


Logical variable for fuel 
Calculation 


Make fuel calculation 


Do not make fuel calcula- 
tion 





Variable 
PARM (3) 

ENGN (3, IENGN) 
EHP (I) 

DEL (3) 

DELMIN (3) 

VK (I) 

IPRNT 

DISPL 


IN 
MAX 


XJ 
IK 
IL 


KM 
Kl 


WEIGT 


Te 


DISPL 
I 


IJ 


VARIABLES USED IN H2OJT 
NOT IN COMMON 


Description 
Working system parameters 
Contains name of engine IENGN for output 
Effective horsepower at condition I, horsepower 
Initial step sizes of system parameters 
Minimum step sizes of system parameters 
Craft speed at condition I, knot 
Print data set reference number 
Craft displacement, 


long ton 


Estimate of minimum number of gas turbines 
required to power craft 


Estimate of maximum number of gas turbines 
needed to power craft 

Working variable for number of gas turbines 
Working variable for craft condition I 
Working variable for IN 


Minimum number of pumps to be used for XJ 
gas turbines 


Maximum number of pumps to be used for XJ 
gas turbines | 
Propulsion system weight, long ton 


Working variable for prime mover character- 
istic L | 


Craft displacement, long ton 


Index for do loops 


Working variable for check on power adequacy 


Index for do loop 
-84- 





Variable 


ISAVE 


N 


XMIN 


bescription 


Working variable for MAX 


Working variable for output of. prime mover 
type 


Factor for minimum acceptable prime mover 
SHP operation 


Working variable for prime mover character- 
istic L 
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Variable 


APOD (I) 


ASTRT (I) 


ASPRY (I) 


211) 


TSUM 


DSIG 
TOTAL 
SIGMA 


SUMZ 


SUMX 


SUM 


250М 


XSUM 


SPD 


JDRAG 


LIST OF VARIABLES USED IN FCT 


NOT IN COMMON 
Description 

Difference between present and previous 
nacelle drag calculations at condition I, 
pound force 
Difference between present and previous strut 
draq calculations at condition I, pound 
force 
Difference between present and previous spray 
drag calculations at condition I, vound 
force 
Margin factor for thrust at condition I 


Workina variable for propulsion system 
weight, pound 


Workina variable for SIGMA 
Total head less vapor head, foot 
Head above vapor head, foot 


Working variable for moment in vertical plane, 
due to weights, foot pound force 


Working variable for moment in longitudinal 
plane due to weights, foot pound force 


Working variable for system weights, pound 


Total moment in vertical plane due to weights, 
foot pound force 


Total moment in lonaitudinal plane due to 
weights, foot pound force 


Working variable for jet velocity calcula- 
tion, foot per second 


Working variable for NUMB 
Indicator of acceptable drag accuracy 


New drag calculation is within 5% of previous 
Ре 





Variable 


ANGLE 


COLF 


KOUNT 
T 


PARM (K) 


DDRAG (I) 


BSPRY (I) 


DSTRT (1) 


DPOD (I) 


Description 
calculation 
New drag calculation is greater than 5$ of 
previous calculation, redesign strut and 
nacelle 
Optimum nozzle depression angle, radian 
Factor for nacelle drag calculation 
Strut must be resized due to cavitation 
Use nacelle design 
Working variable for craft condition I 


Craft condition 


Working variable for system parameters 


VJVO 
VIVO 


Previous acceptable total drag calculation at 
condition I, pound force 


Previous acceptable spray drag calculation at 
condition I, pound force 


Previous acceptable strut drag calculation at 
condition I, pound ferce 


Previous acceptable nacelle drag calculation 
at condition I, pound force 
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Variable 


ZK 


SPO 


РАР 


SIGTV 


JNUMB 


ЕГЕХ 


XD 


DIDMX 
01 

DI 

CP IN 


QIN 


QC 


QAUX 


KDEX 


VR 


LIST OF VARIABLES USED IN NACEL 


NOT IN COTTON 
Description: 


Percentage of auxiliary inlet area permitting 
flow 


Static pressure at inlet, pound per square 
foot 


Vapor pressure, pound per square foot 


Incipient turnina vane cavitation number, 
referenced to diffuser exit velocity 


Working variable for NUMB 
Craft condition, also index for do loops 
Peak external pressure coefficient 


Workina index for do loop for ELEXT/DM 
values | 


Working index for do loop for VIVO values 


Length of forebody to maximum diameter ratio, 
ELEXT/DM 


Maximum permissible diameter ratio 

Flow rate ver strut, cubic foot per second 
Nacelle inlet diameter, foot 

Peak internal pressure coefficient 


Take-off flow rate per nacelle through inlet, 
cubic foot per second 


Take-off flow rate per nacelle, total, cubic 
foot per second 


Auxiliary flow rate required to avoid 
cavitation, cubic foot per second 


Counter for iterations on diffuser length 
Velocity ratio at take-off condition, based 


on flow through inlet 
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Variable 


ү12 


PRL2 


ELI 


SPI 


VIAUX 


PRAUX 


PTAUX 


DYP 


РС 


ODIF 


ЕНІ 


PHS 


D2 
01 
ELMAX 
ELMIN 
II 
EL 


DEL 


Description 


Inlet velocity at take-off condition, 
based on flow throuch inlet 


Pressure recovery coefficient of lip at 
take-off 


Inlet stagnation pressure immediately aft 
of lip at take-off, pound per square foot 


Inlet static pressure immediately aft of lip 
at take-off, pound per square foot 


Inlet velocity in auxiliary inlet, foot per 
second 


Pressure recovery of auxiliary inlet 


Staanation pressure inside auxiliary inlet, 
pound per square foot 


Net dynamic pressure immediately aft of lip, 
pound per square foot 


Average inlet staanation pressure of combined 
flow, pound per sauare foot 


Working variable for diffuser flow rate, 
cubic foot per second 


Equivalent angle of forebody, radian 

Sine of equivalent angle of forebody 

Lenath of auxiliary inlet, foot 

Diffuser exit diameter, foot 

Inlet diameter, foot 

Maximum permissible length of diffuser, foot 
Minimum permissible length of diffuser, foot 
Working variable for craft condition I 
Working variable for diffuser length, foot 


Workina variable for change in diffuser lenath, 
foot 80 





Variable 


ELD 
ELL 


ELFAC 


DDM 


XKT 


DL 


CDRG 
ANGL 
EDIF 


POW 
POWI 


EM 
AEXN 
REND 


DDIF 
PLOSS 


VAOUT 
SQUAR 


VCRIT 


Description 
Working variable for nacelle length, foot 
Working variable for nacelle length, foot 


Working variable for excess in nacelle length 
due to diffuser, foot 


Average diffuser diameter, foot 

Form loss coefficient of diffuser 
Reynolds number, based on nacelle length 
Reynolds number, based on inlet diameter 


Ratio of maximum external diameter to nacelle 
length 


Computed drag coefficient 
Equivalent half angle of diffuser, degree 
Diffuser expansion factor 


Power loss due to drag and duct loss of 
diffuser, horsepower 


Previous power loss calculation for diffuser 
length, horsepower 


Factor in wetted surface calculation 
Wetted surface area, square foot 
Reynolds number, based on inlet diameter 


Total pressure loss in diffuser, pound per 
square foot 


Total pressure loss in nacelle, pound per 
square foot 


Average exit velocity, foot per second 
Factor in critical velocity calculation 


Critical velocity in strut elbow, foot per 


second 
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Variable 


VMAX 


RENL 
NL(2) 
ML (2) 
KL (2) 
JL (2) 
IL (2) 
VPT (6) 


XDT (10) 


PRLT (6) 
SIGI (I) 


PTO(I) 


XDTT (J) 


CDUMX (K) 
CDUMY (L) 


CD(I) 
DIDMT (J) 


УВТТ (9) 


DLIP (I) 


QO (I) 


VELR (I) 


Description 


Maximum velocity at nacelle exit, foot per 
second 


Reynolds number, based on nacelle length 
Interpolation parameters 

Interpolation parameters 

Interpolation parameters 

Interpolation parameters 

Interpolation parameters 

Data array of velocity ratios 


Tabulated forebody length to inlet diameter 
ratios 


Tabulated lip pressure recovery coefficients 
Free stream cavitation index at condition I 


Stagnation pressure at craft condition I, 
pound ver square foot 


Tabulated ELEXT/DM ratios for trim angle J 


Dummy array of peak pressure coefficients for 
VIVO K 


Dummy array of velocity ratios for ELEXT/DM 
L 


Drag coefficient at craft condition I 
Tabulated DI/DM values for ELEXT/DM J 


Dummy array of velocity ratios for angle of 
attack (J), internal 


Lip loss coefficient at craft condition I 


Free stream dynamic pressure at craft condition 
I, pound per square foot 


Inlet velocitv ratio at craft condition I 
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Variable Description 


VRTEX (I) Same as VRTT(I) but external 

VRMAX (1) Maximum permissible velocity ratio at craft 
condition I 

FRICT Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 


dioe 





Variable 
SHAPE (3,ITYPE) 
THETA(J) 


XLOSS (J) 


ROA (10) 
RE 


ERICT 


CORR 


KOUNT 


IELB 


FACTR 
RIN 

ROUT 
RATIO 


XN 


VARIABLES USED IN ELBOW 
NOT IN COMMON 


Description 
Contains name of shape ITYPE 
Data array of elbow angles 


Data array of elbow loss coefficients with 
thin, circular arc turning vanes 


Outside radius of splitters, foot 

Reynolds number of duct 

Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 

Function statement,.calculates head loss 
correction factor for different Reynolds 


numbers | 


Maximum Reynolds number permitted for 
splitter loss equation 


Working variable for craft condition I 


Index indicating which elbow is being 
designed/evaluated 


Strut 

Hull 

Pump 

Divergence 

Factor used in splitter loss calculation 
Inside radius of bend, foot 

Outside radius of Беко, foot 

Desired radius ratio, 4.3 


Number of subdivided elbows required to 


achieve BE cm 





Variable 


N 


М1 


SUM 


XCORR 


ЖКТ 


DIAM 


VOLV 


VOL 


AREAL 


IJ 


IK 


HEAD 


Description 
Number of subdivided elbows used 
Number of splitters corresponding to N 


Working variable of sum of head times 
subdivided elbow area, foot cubed 


Inside radius of subdivided elbow, foot 
Average velocity, foot per second 
Heiaht of subdivided elbow, foot 


Equivalent cross-sectional area of subdivided 
elbow based on HGT, square foot 


Equivalent radius of subdivided elbow, foot 


Ratio of head loss correction factors for 
differing Reynolds numbers 


Head loss coefficient of subdivided elbow 
Equivalent diameter, foot 
Volume of splitters, cubic foot 


Volume of splitters and elbow structure, 
cubic foot 


Duct area ver elbow, square foot 
Working variable for shape determination 
Working variable for shape determination 


Working variable for number of subdivided 
elbows N 


Ratio of inside radius to outside radius 
of subdivided elbows 


Average head loss in elbow, foot 
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Variable 


ШНЕТ2 


ЕХРАМ 
ARATO 


VL 


SIGMA 
CM 
DEIN 
WIDE 
DEOUT 
DEAVE 
SIRT 
ALONG 
STAN 
THETA 
E COEF 
FORML 
VELIN 


VLOUT 


CDS 
CDSP 


RE 


LIST OF VARIABLES USED IN STRUT 


NOT IN COMMON 


Description 


Array of data of equivalent angle of diffuser, 
degree 


Array of data of expansion coefficient 
Area ratio of strut diffuser 


Local maximum velocitv, external, foot per 
second 


Local cavitation number, external 

Mean chord, foot 

Strut inlet @egqgulvalent diameter, foot 
Width of duct at strut exit, foot 

Strut exit equivalent diameter, foot 
Average equivalent diameter, foot 
Vertical strut length, foot 

Actual strut length, foot 

Arctangent of equivalent diffuser angle 
Equivalent angle of diffuser, 20, degree 
Diffuser expansion factor 

Diffuser expansion loss coefficient 
Average inlet velocity, foot per second 
Average exit velocity, foot per second 
Strut Reynolds number, based on mean chord 
Strut drag coefficient 

Spray drag coefficient 


Duct Peynolds nurber, based on inlet velocity 
E 





Variable 


PIPEL 


TOTAL 


HEAD 


KOUNT 


HGT 


CGWS 


Description 


Duct friction lessi coefficient 

Total loss coefficient 

Head loss of diffuser, foot 

Working variable for craft condition I 
Elevation of strut, foot 


Vertical center of aravity of duct, foot 


296 





LIST OF VARIABLES USED IN JUNCT 
NOT IN COMMON 


Variable Description 

ALPHA Assumed junction angle of 0°, degree 

DIAM Duct diameter, foot 

FRCTL Friction loss coefficient for athwartships 
length 

KOUNT Working variable for craft condition I 

AJCT Length of fore and aft ducting to pump inlet, 
Tot: 

FRCTJ Friction loss coefficient for junction 

XLAMD Working variable for AMIXL calculation 

FLONG Length of athwartships ducting, +7 

AJCTL Total loss coefficient of junction 

AMIXL Mixing loss coefficient of junction 

V Average velocity, foot per second 

XPUMP Number of gas turbines + 1 

RE Reynolds number of duct 

BETA (J) Data array of junction angles, degree 

ALAMD (J) Data array of mixing loss coefficient 


corresponding to BETA(J) 


БЕТСТ Function statement, calculates Moodv friction 
factor for smooth pipe, based on Reynolds 
number 
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Variable 


KOUNT 
V 
DIAM 
RE 


APIPE 
ЕСТІ. 


E PUMP 


XKT 


XLONG 


NT OT 


LIST OF VARIABLES USED IN PIPE 


NOT IN COMMON 
Description 
Working variable for craft condition I 
Average velocity, foot per second 
Duct diameter, foot 
Reynolds number of duct 


Length of fore and aft pipe to pump inlet, 
foot 


Friction loss coefficient for athwartships 
lenath 


Number of gas turbines 


Friction loss coefficient of fore and aft 
length 


Length of athwartships pipe, foot 
Function statement, calculates Moody friction 


factor for smooth pipe, based on Reynolds 
number 


ов 





Variable 


ERICT 


KOUNT 
RE 


DIVL 


DWGT1 
DWGT2 
DWGT3 
WWGT 1 
WWGT2 
WWGT 3 


CGWX1 


CGWX 2 


CGWX3 


DIAM 


FRCTL 


DIVLC 
XPUMP 
ANGLE 
HEADL 


FLONG 


LIST OF VARIABLES USED IN DIVRG 


NOT IN COMMON 
Descrintion: 
Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 
Working variable for craft condition I 


Reynolds number of duct 


Divergence loss coefficient, not including 
friction 


Duct weight of divergence angle, pound 
Duct weight of divergence length, pound 
Duct weight of pump inlet anale, pound 
Water weight of divergence angle, bound 
Water weight of divergence length, pound 
Water weight of pump inlet angle, pound 


Longitudinal center of gravity of divergence 
angle, from transom, foot 


Longitudinal center of gravity of divergence 
length, from transom, foot 


Longitudinal center of gravity of pump inlet 
angle, from transom, foot 


Duct diameter, foot 


Friction loss coefficient for athwartships 
length 


Total divergence loss coefficient 

Number of gas turbines + 1 

Divergence angle, also pump inlet angle, radian 
Total divergence head loss, foot 


Athwartships length, foot 
EU 





Variable 


XLONG 

V 

ADIV 
THETA (J) 


COEF (J) 


Description 
Fore and aft Wengen, foot 


Average velocity, foot per second 
Divergence length, ADIV = XLONG/COS (ANGLE), foot 
Data array of divergence angles, degree 


Data array of divergence loss coefficients 
without friction, corresponding to THETA(J) 


= 1002 





Variable 


KOUNT 
XLPS 


XFAC 


XPUMP 
ANOZ 
XLNOZ 
DT 

DJ 


XCORR 


АББА1 
AJETI 


QQ 


ЕЕГ СТ 


LIST OF VARIABLES USED IN NOZZL 


NOT IN COMMON 
Description 
working variable for craft condition I 
Dummy nozzle length, foot 


Check on whether nozzle exits through bottom 
or stern 


Number of pumps 

Optimum nozzle depression angle, radian 
Nozzle length, foot 

Nozzle throat diameter, foot 

Nozzle jet diameter, foot 

Nozzle head loss factor 


Reynolds number, based on average diameter 
and velocity 


Throat pipe area, square foot 

Jet pipe area, square foot 

Flow rate per nozzle, cubic foot ver second 
Function statement, calculates Moody friction 


factor for smooth pipe, based on Reynolds 
number 
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LIST OF VARIABLES USED IN PUMP 
NOT IN COMMON 


Variable Descrivtion 


K Cruise condition indicator 

J Take-off condition indicator 

I Index for do loops 

XP UMP Number of pumps 

THOMI Lower limit of Thoma's criterion for sinale 


inducer axial pump 


DRAT Hub to tip diameter ratio 

XXLP Factor for pump length 

CW Weight coefficient 

QX Ratio of cruise flow rate to take-off flow 
rate 

CA Factor for RPM calculation 

CB Factor for RPM calculation 

HX Patio of design pump head to off design head 

RX Ratio of design to off design pump RPM 

PHRAT Ratio of off desian to design flow 
coefficients 

XNGT Number of das turbines 

ETAPP Product of pump and gearbox efficiencies 

HP ہج‎ Inducer head for one and two stage axial 


pump designs, foot 


THOMS Phona' swavitation eriterion for inlet to 
axial stage 


HHP Axial stage pump head, foot 


BETA2 Exit жес» radian 
ЖАР 





Variable 
XNNS 

BD 

ЕС 

ДА 


IMPL 


N 
РС(К,Ј) 
РСА(К,Ј) 
ШЕС(К, 7) 


XRPM(M) 
XD1 (M) 
XD2 (M) 
RPK (M) 
XPUP (M) 
YLP (M) 
XERAT (M) 


APUP (NSTG) 


WRAT (NSTG) 


Description 
Non-dimensional specific speed 
Impeller exit width ratio 
Paetor met lew coefficient calculation 
Factor in flow coefficient calculation 


Maximum number of impellers permitted for 
centri fugal purp, = 10 


Workina variable for number of centrifugal 
oump impellers 


Working variable for NSTG 

Inducer head curve coefficients 

Inducer plus axial stage head curve coefficients 
Centrifugal pump head curve coefficients 


Working variable for off design RPM of 
centrifugal pump with M impellers, RPM 


Working variable for inlet tip diameter, DIS, 
foot 


Working variable for exit tip diameter, D2S, 
toot 


Working variable for desian RPM of centrifugal 
pump with M impellers, square foot 


Working variable for centrifugal pump area 
with M imvellers, foot squared 


Working variable for centrifugal pump length 
with M impellers, foot 


Working variable for gear ratio required for 
centrifugal pump with M impellers 


Inlet area of pump NSTG, square foot 


Weight ratio of pumo NSTG, including pump dry 
and wet weight, gearbox and fuel 


=) з= 





Variable Description 


WD (M) Working variable for XWD(NSTG) for 
centrifugal pump with M impellers 


WW (M) Working variable for XWW(NSTG) for centri- 
fugal pump with M impellers 


WWG(M) Workina variable for XWG(NSTG) for centri- 
fugal pump with M impellers 
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Variables 
CES 
CA 


CD 


XN 


TI 


IJK 


XJ 


ОПРР 
SHNG 
M 


WT (ITIME) 


DIS(ITIME) 


VJJ (ITIME) 


LIST OF VARIABLES USED IN FUEL 


NOT IN COMMON 
Description 
Constant for SFC «calculation 
Cruise condition drag to lift ratio 


1 + total system head loss coefficient, 
based on jet velocity 


Number of intervals endurance is divided into 
+ ] 


Number of intervals endurance is divided into 


Time to cover one ranae interval at constant 
VOU), hour 


Index for shifting SFC curves 


Factor to convert SFC if SHP is less than 70% 
of design SHP 


Working index for N 

Head required, foot 

Total thrust required, horsepower 

Total thrust required per engine, horsepower 
Working index for N 


Weight of fuel used in time increment ITIME, 
pound 


Displacement at time increment ITIME, pound 


Jet velocity at time increment ITIME, foot 
per second 


EIE 





CISIT OP VARIABLES USED IN PTIRN 
NOT IN COMMON 


Variable Description 

PSI (N) Current basepoint coordinate of parameter N 
THETA (N) Previous basepoint coordinate of parameter N 
PHI (N) Present exploratory point coordinate of 


parameter N 


DEL (N) Current step size of parameter N 

DELMIN (N) Minimum step size of parameter N 

DIR(N) Last successful direction of parameter N 

SAVE (N) Working variable for PHI (N) 

5 Working variable for function value 

SPHI Working variable for function value at PHI 
coordinates 

SPSI Current best function value at PSI coor- 
dinates 

NUMB Counter for minimum step size check 

RHO Step size change factor 

ICALL Indicator of current point move 

K Index for do loop 

I Index for do loop 

N Number of parameters of search 

SIGN (N) Directed step size of parameter N 
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LIST OF VARIABLES USED IN OUTPUT 


Variable 
IPRNT 
J 


SUM 


DRATO 
IK 


XNGT 


LL 


KLL 


NIMP 

TFM 

WT RAT 

WRATF 

ENGN (IENGN) 
VJRAT (I) 
VIRAT ( I) 


HEADL (I,ICOMP) 


CONDS (2,1) 


NOT IN COMMON 
Description 
Print data set reference number 
Working variable for condition I 


Total duct head loss, excluding elevation, 
foot 


Index for do loop 

Strut diffuser area ratio 

Index for number of elbows in system 

Number of gas turbines 

Index for do loop 

Index for implied do loop in output statement 


Working variable to point to correct head 
loss for output 


Working variable to point to correct format 
statement for head loss output 


Index for implied do loop in output statement 
Number of impellers in pump 

Strut thickness at flying waterline, foot 
Total propulsion system weight ratio 
Propulsion system weight ratio, excluding fuel 
Contains name of engine IENGN 

Jet velocity ratio at craft condition I 

Inlet velocity ratio at craft condition I 


Head loss of component ICOMP at craft con- 
dition I 


Label for craft condition I 
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Variable Description 


ELBWS (2,1K) Label for elbow IK 
РС(Т) Propulsive coefficient at craft condition I 
LABEL (5,M) Labels for output 
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